Isocurvature, non-gaussianity and the curvaton model 
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Recent analyses of the statistical distribution of the temperature anisotropies in the CMB do 
not rule out the possibility that there is a large non-gaussian contribution to the primordial power 
spectrum. This fact motivates the re-analysis of the curvaton scenario, paying special attention 
to the compatibility of large non-gaussianity of the local type with the current detection limits 
on the isocurvature amplitude in the CMB. We find that if the curvaton mechanism generates a 
primordial power spectrum with an important non-gaussian component, any residual isocurvature 
imprint originated by the curvaton, would have an amplitude too big to be compatible with the 
current bounds. This implies that the isocurvature mode should be equal to zero in this scenario 
and we explore the consequences of this inference. In order to prevent the generation of a such a 
signal, the CDM must be created at a late stage, after the curvaton decays completely. This is used 
to constrain the nature of the CDM, arriving at a general relation between the temperature of the 
universe at CDM creation and the scale of inflation. It is possible to find an absolute maximum for 
the temperature at CDM creation, which is dependent on the particular infiationary potential. For 
a quadratic potential, we find Tcdm < 1.7 x lO^GeV. 



I. INTRODUCTION 

An analysis performed by Yadav and Wandelt [Tj, 
claims the detection of a significant amount of primor- 
dial non-gaussianity of the local type in the WMAP 3- 
year data. More recently, the WMAP team [2| published 
their latest data release in which the presence of a non- 
gaussian component in the primordial power spectrum is 
analyzed. By using a more sophisticated and conserva- 
tive mask, the bounds found for the non-gaussian con- 
tribution are again reconciled with no presence at all, 
however the region of a high non-gaussian contribution 
has not been ruled out. 

The clear detection of such a signal would automati- 
cally rule out the simple single-field inflationary model, 
since one of its predictions is a gaussian spectrum for the 
density fluctuations. We believe that there is not strong 
evidence against the existence of such a signal and there- 
fore, it is interesting to investigate more sophisticated 
models that could explain the origin of the primordial 
anisotropies along with primordial non-gaussianity at the 
levels detected. 

One suitable candidate, is the curvaton model [SlU- 
It does predict an almost flat adiabatic spectrum for the 
density fluctuations along with a possibly large, local, 
non-gaussian signal and, also possibly, a wholly corre- 
lated isocurvature signal which is inside the allowed ex- 
perimental range up to now. 

While we wait for the confirmation or rejection of the 
non-gaussian signal, we study the new constraints im- 
posed in a general curvaton scenario in which the non- 
gaussian component contributes to the total power spec- 
trum at the level of about 0.1%. In particular, we study 
the compatibility of an eventual detection with the cur- 
rent bounds on primordial isocurvature and see what con- 
sequences this may bring to the cosmological model. 

This analysis differs that that carried out in [2] in the 
sense that they only examine the curvaton as a possible 
source of isocurvature whereas we study the possibility 



that this mechanism generate a large non-gaussian signal 
and maybe, some residual isocurvature. The approaching 
angle of this work is more similar to what is presented in 
[5] although different assumptions and new data, allow 
us to draw independent results. 

In section |ll] we outline the curvaton model to arrive 
at an expression for the predicted non-gaussianity and 
isocurvature signals in terms of the factor r, the ratio of 
the energy density of the curvaton to the total energy 
density at the decay of the curvaton. In section |III| we 
present the results from several analyses on the depar- 
tures from gaussianity and adiabaticity and find which 
boundaries are imposed on r by each one of them. In 
section |IV| we combine those constraints to find interest- 
ing implications for the decoupling temperature of a cold 
dark matter candidate. Finally, we present our conclu- 
sions section IVl 



II. THE CURVATON MODEL 

The curvaton infiationary model, was proposed as an 
alternative to the usual mechanism for seeding the pri- 
mordial curvature perturbations [3l |4]. Rather than 
evolving from quantum fiuctuations in the infiaton field, 
curvature perturbations arise due to the presence of an 
aditional scalar field, the curvaton (a). This field is prac- 
tically massless and its energy density, pa, is highly sub- 
dominant at very early times so it is a spectator during 
infiation. However, after infiation ends and the infiaton 
decays completely into radiation, pa- becomes more and 
more relevant and the fiuctuations in the curvaton field 
emerge as the predominant seed of the observed struc- 
ture. After that, the Hubble factor decreases to a value 
close to the curvaton mass, and a starts oscillating. It 
finally decays and the regular Hot Big Bang evolution 
proceeds. 

There are three general conditions that must be ful- 
filled in order for the mechanism to work as described 



above, namely: 

1 . The curvature perturbation generated by quantum 
fluctuations in the inflaton field, C^ , must be negli- 
gible compared to the total curvature perturbation. 

2. The curvaton cannot trigger a second period of in- 
flationary expansion, therefore, its energy density 
must be much smaller than the resulting radiation 
energy density right after inflation. 

3. The curvaton is practically massless during infla- 
tion and thus the value of the field is fixed at ct* . 

(^0 is the curvature perturbation in the uniform density 
gauge, and it is related to the curvature perturbation in 
any other gauge, "01 by: 



which is valid for any generic C,. 



(1) 



A. Details of the model 

We now describe the basic dynamics of the model 
adopting a parameterization close to that used in the 
simplest curvaton model |6, . We summarize here some of 
their results that we will be using for our work and extend 
their analysis by using a general form for the inflationary 
potential and by including the generation of non gaussian 
perturbations. Weusetheset: {T^((/)), to, rcr,o-*}, the po- 
tential for the inflaton fleld and the mass, the decay rate 
and the initial value of the curvaton field respectively. 
Later on, we add the non-gaussian contribution and the 
isocurvature fraction parameters, which will be defined 
in the next sections. 

The total potential during inflation is: 
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which is a good approximation for low values of the cur- 
vaton field. 

We adopt the "curvaton hypothesis" , where the infla- 
ton curvature perturbation is taken to be less than 1% of 
the observed value [7]. Using the COBE normalization 
at the pivot scale, we can set an upper bound for the 
power spectrum of the inflaton: 
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after inflation, condition 2 implies: 

m^cr^ ^ v{(t)) ^ y ~ V{<i)) 



(5) 



As long as the inflationary potential is suitable for a Slow 
Roll hypothesis, we get the upper bound for the value of 
the potential 50 e-folds before inflation ended: 
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H50 < 1.1 X lO'-'GeV 
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Since we do not assume any particular shape for the po- 
tential, we arrive at these bounds by only imposing the 
slow roll condition for the factor e < 1. 

Condition 3 implies that the power spectrum for the 
density fluctuations of the curvaton is given by: 
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where iJ* is the value of the Hubble factor at horizon 
exit. It is shown in [3] that the power spectrum for the 
density contrast for the curvaton field is: 
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where Uexit is the value of the field at horizon exit. Since 
the field is massless during infiation, its value does not 
evolve with time and we fix CTexit — f* for all scales. This 
perturbation spectrum is transferred to the total density 
perturbation by the effects of the non adiabatic pressure: 
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that arises due to the presence of two different, non in- 
teracting, fluids. It is shown that the final value of the 
curvature perturbation inside the sudden decay approxi- 
mation is [5]: 



C^<a 



(10) 



where r is practically the ratio of the curvaton energy 
density to the total energy density at the decay of the 
curvaton, when H ^ T„. In the case where the curvaton 
dominates completely before its decay, r = 1. Using (111), 
([9]) and ( 10 1 , we arrive at the predicted power spectrum 
for the curvature perturbation: 
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which is one of the experimental observables. 

One can write the value of r as a function of the pa- 
rameters of the model. We start with: 
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where the sub and super indices "end" refer to the value 
of the variables at the end of infiation and am is the scale 



factor at the time when the curvaton starts oscillating, 
when H r^ va. Assuming that the decay occurs during a 
phase of radiation domination (r < 1), 
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Plugging in the values of the energy densities at the end 
of inflation, we arrive at: 



(13) 
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Note that the expresion above is independent of the as- 
sumed potential for the inflaton. 

In order to prevent a second period of inflation, the 
curvaton must provide a negligible contribution to the 
energy density, at least until it starts oscillating, at iJ ~ 
m. This implies: 



If we identify the first term of ( 19 ) as <I>l, then, compar 



ing to (16 1 we get [9j: 
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Note that since we have not taken into account the intrin- 
sic, second order perturbations in the Bardeen potential, 
this expression is only valid for high values of /nl which 
dominate over the intrinsic non-gaussianity. This approx- 
imation is good enough for our purposes since, motivated 
by the claim in [T] , we assume study the consequences of 
a local /nl '^-^ 100 (which makes the non-gaussian sig- 
nal contribute with 0.1% to the total curvature pertur- 
bation). A more precise definition of the non-gaussian 
contribution can be found in |11) . 

The assumed value for /nl implies a value of r = 
0.0125 which also validates the hypothesis assumed in 
the previous section that the curvaton decays during ra- 
diation domination. 
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On the other hand, the decay must occur after the oscil- 
lation of the curvaton: 



Trr < m 



(15) 



The two factors of equation ( 13 1 have opposite effects 



that prevent us from making any a priori estimation of 
what the value of r might be. Fortunately, r is related 
to two other observables apart from the power spectrum: 
the amount of non-gaussianity, /nl and the isocurvature 
amplitude, fiso- This will help us put tighter bounds on 
its value. 



B. Non-gaussianity in the curvaton scenario 

The magnitude of the non-gaussian contribution, is 
parametrized as pUj : 



C. Isocurvature in the curvaton scenario 

We focus on the generation of a CDM-isocurvature 
(GDI) amplitude inside the curvaton scenario. 
The GDI mode is defined as: 



^cdni — 



Sp, 



cdm 



Pcdr, 



3 Spj 

4 p^ 



(21) 



This quantity is gauge invariant and, in terms of the in- 
dividual curvature perturbations (in the flat gauge): 
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We assume that all the species in the universe are radi- 
ation when the GDM is created, therefore, we can write 
C^ = (■. As in ref. [9], we take the epoch of creation of the 
dark matter particle as the time after which its number 
density is conserved. 

There are several parameter configurations that could 
potentially lead to the generation a GDI mode, depending 
on when the GDM is created: 



* = *L + /NL*i 



(16) 



where <& is the Bardeen potential and is related to the 
curvature perturbation as: 
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for super-horizon modes during the matter domination 
epoch. We can expand pa- up to quadratic order: 



1. The CDM IS created after the curvaton decays completely 

In this scenario no isocurvature is generated. The 
species from which the GDM is hypothetically created, 
should have thermalized with the rest of the universe, 
carrying the same curvature perturbation, (^. This is the 
perturbation that will be inherited by the GDM fiuid so 
that Ccdm = C- Thus, (22 1 is zero. 



Pa = mla^ 



and plug it into (17 1 
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2. The CDM is created before the curvaton decays and 
while f <^ 1 




(19) 



In this scenario, Ccdm = and according to ( 22 1 : 
'^cdm = ^3C 



(23) 



that is, the isocurvature amplitude is totally anti- 
correlated to the adiabatic mode, and three times as big. 
This possibility has been ruled out by the analysis of sev- 
eral datasets performed by various groups [T^IOirH] . 



3. The CDM ts created while the curvaton decays 

This scenario is the most model dependent, and we 
distinguish two possibilities: i) the CDM is created while 
the curvaton decays and pa- is comparable to pr- It is pre- 
dicted that some GDI will be generated, but the precise 
amount depends on the particular creation mechanism. 
This has been left for further study, ii) the curvaton 
decays into the CDM. Then, the CDM curvature per- 
turbation will be Ccdm — Ccti E^nd using (10 1, we arrive 
at: 



^cdm — 3 ( ) C 



(24) 



Which leads to an isocurvature mode that is correlated 
to the adiabatic one, has the same scale dependence, and 
has an amplitude that depends on the ratio of curvaton 
to radiation energy density. The isocurvature amplitude 
is small if the curvaton is close to dominating the energy 
density of the universe when it decays. If the curvaton 
decays while it is still sub-dominant, the resulting en- 
tropy mode could be large. 

In section [IlB I we mentioned the value of r correspond- 
ing to a 0.1% non-gaussian contribution in the primordial 
power spectrum: r = 0.0125. In what follows, we see that 
the current upper bounds on a curvaton-like isocurvature 
signal are in tension with such a value, that indicates 
that the curvaton energy density is not close to domi- 
nation when it decays. This starts pointing towards our 
conclusion that if the curvaton mechanism is responsible 
for a large non-gaussianity , then, the same mechanism 
cannot generate any isocurvature signal whatsoever. We 
see this in more detail in the next section. 



More recently, with a better understanding of the point 
source contamination and a more restrictive KQ75 mask 
[TB], the WMAP team reported an allowed region [5]: 



9 < /nl < 111 (at 95% c.l.) 



(26) 



from the analysis of the bispectrum and using the same 
estimator as Yadav and Wandelt. The maximum mul- 
tipole number used in [J is fmax = 700. The WMAP 
team also reports a second result from the analysis of the 
Minkowski functionals following the method described in 
[T7] . With this method, they find a new result which is 
in slight tension with the previous one: 



178 < /nl < 64 (at 95% c.l.) 



(27) 



Recently, Hikage et. al also studied the non-gaussianity 
in the WMAP 3rd-year data with the Minkowski func- 
tionals 118 to find: 



70 < /nl < 91 (at 95% c.l. 



(28) 



which seems in agreement with the WMAP team analysis 
with the Minkowski functionals. 

We show all these different results to highlight two 
points: one is the fact that there is al lot left to under- 
stand about primordial non-gaussianity and the statisti- 
cal methods to detect it. The other is that none of the 
results quoted above rule out a possibly large (of order 
0.1% ) non gaussian component and that is the reason 
why we believe it is worth exploring the curvaton model 
as the responsible for a large /nl- 



Note that a part of the two sigma allowed region in ( 26 1 



cannot be explained in terms of the curvaton dynamics 
because r has a natural upper bound of 1 . This value of 
r corresponds to /™l" = —1.25 (taking into account the 
second order corrections to /nl described in [11_). 



B. Isocurvature 



III. DATA ANALYSIS 

In this section we describe the existent experimental 
constraints on /nl and | "?*" | . The most restrictive re- 
sults to the date are those reported by the WMAP team 
in their 5- year report, which we quote along with signif- 
icant analyses performed by other groups. 



A. Non-gaussianity 

In [T], Yadav and Wandelt analysed the two main 
science channels of the WMAP 3rd-year data up to 
^max = 750. Using the estimator described in [15; and 
the KpO foreground mask, they claimed a detection of a 
non-gaussian, local, component of: 



26.91 < /nl < 146.71 (at 95% c.l. 



(25) 



With two more years of data gathering the WMAP 
team achieved stringent limits on the parameters of the 
basic cosmological model and some extensions such as 
departure from adiabaticity. They use the parameter a 
to measure the amplitude of the isocurvature mode in 
the scenario in which the curvaton decays into dark mat- 



ter generating a mode with amplitude described by ( 24 1 
The parameter a is defined as [19] : 



J ISO 



l + /i 



2 
iso 



(29) 



where /iso is the ratio the isocurvature and adiabatic am- 
plitudes at the pivot scale. In the notation defined in the 
sections above, /iso = | 7?"* | . This parametrization has 
the advantage that a is naturally bounded to take a value 
e [0,1]. Some controversy exists about whether or not 
some isocurvature models are artificially favored by this 



choice [m [20] but our model is already selected, so this 
should not be an issue for the analysis. 

The WMAP team finds the following upper bound for 
a inside the frame of the curvaton model (totally cor- 
related ^, same scale dependence of the adiabatic and 
isocurvature amplitudes): 



< 0.0037 (at 95% c.l.) 



(30) 



using their 5-year CMB data plus Type la supernovae 
[H] and Baryon Acoustic Oscillation data [22]. This 
bound improves previous results obtained with for the 
3rd-year WMAP data release plus other data, including 
large scale structure [T^ by a factor 2.5. However, the 
tightest bound to the date is found when the Lyman-a 
forest data is included, a < 0.0015 (at 95% c.l.) ^. 

In any case, the difference in the derived value of the 
ratio of energy densities is minimal. Taking into account 
(24 1, we see: 
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Thus, we finally arrive at the experimental result that if 
the curvaton decays into dark matter, it must do so when 
its energy density is just about to become the dominant 
component in the universe, i.e., the value of r must be 
bounded by: 



0.98 < r < 1 



(32) 



(using the WMAP bound on a) . This bound corresponds 
to a non-gaussianity contribution of: —1.21 > /nl > 
— 1.25 which is clearly inside the 2 a region derived by 
f2]. However, it is one sigma away from the central value 
quoted in Table 5 of [2J. Although this fact is not signif- 
icant from a statistical point of view, it suggests the in- 
triguing possibility that the non-gaussianity values could 
be in strong confiict with the value of /nl derived from 
the curvaton-like isocurvature bounds. 

In what follows, we focus in this potential scenario, and 
study the consequences of an eventual detection of a high 
non-gaussian component inside the curvaton framework. 



IV. IMPLICATIONS FOR CDM 

As described in the previous section, the detection of 
a large non-gaussianity would restrict the allowed bound 
for the ratio r to a region incompatible with the range 
suggested by current limits on the isocurvature ampli- 
tude. 



^ In [2], they use the "anti-correlated" result, but this is due to 
the sign convention used for the definition of the correlation am- 
plitude. 



Thus, if one accepts that the curvaton is the model re- 
sponsible for the generation of the primordial non gaus- 
sianity, one must accept as well the fact that the curvaton 
mechanism proceeds according to subsection |IIC 1| and 
no isocurvature imprint is left in the CMB. This leads to 
two important consequences: 

1. The CDM is not the direct decay product of the 
curvaton. 

2. The CDM must be created after the curvaton field 
has decayed. 

The last point implies a severe constraint on the possible 
values of the temperature of the universe at the epoch of 
creation of dark matter. For thermal relics, this trans- 
lates directly onto a constraint on their mass. Inversely, 
the detection of a particular dark matter species, would 
give us some hints about the precise dynamics of the cur- 
vaton mechanism. 

Firstly, let us point out that consequence 2 implies: 



Hr 



<r. 



(33) 



where Hdm — 2^c^dm/-^'^PL is the Hubble factor at the 
epoch of DM creation, at a temperature Tcdm- The upper 
bound for Tcdm is then. 



Tcdm < yr.M] 



PL 



(34) 



Thus, a constraint in T^- could have strong implications 

for the particular mechanism of dark matter generation. 

Note that the possibility that the CDM was created 

much before the curvaton decayed, had been ruled out 



already, as we mention in section II C However, this fact 



alone was not enough to impose severe constraints on 
Tcdm- In the first place, one could always appeal to the 
attractive possibility that the curvaton decayed indeed 
into the dark matter- an instance that has been ruled 
out under our assumptions. Secondly, without the ad- 
ditional constraint provided by the assumption of large 
non-gaussianity, it is not possible to derive the tight 
boundaries on the decay epoch of the curvaton that we 
describe below. Therefore, it is this particular combina- 
tion of extensions to the standard inflationary scenario 
that allows us to draw the conclusions that we present as 
our main result. 

In what follows, we go through further considerations 
about the model that help us reduce the dimension of 
the space of parameters and arrive at a simple relation 
between the maximum Tcdm and the scale of inflation. 

We can solve for F^ in (13 1 and write r as a function 
of /nl: 



F. 



/5/nl cr* 

V 4 
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(35) 



Its value depends on m and ct* which are in principle very 
loosely bounded. Below we will show how two additional 
considerations relate the value of two of these parameters 
to each other. 



The value of ct. 



Even though the bound ( 14 1 trims off a large range 



of possible initial values of the curvaton field, it is not 
informative enough for the purpose of studying the limits 
of the mass of a cold dark matter candidate. Luckily, it 
is possible to find a lower bound for a^ if we consider the 
effect of the de Sitter vacuum on the field. It is shown in 
[23] that the value of the variance of an almost massless 
field after a period of exponential expansion approaches 
the limit: 



(^^) = 



SiJ-* 



(36) 



if the expansion has occurred for a long enough time (i.e., 
the number of e-folds before iV ~ 50 is large). It is 
sensible then, to expect that the field at least has moved 
that much away from its minimum, and thus, this value 
can be taken as the lower bound for cr* after inflation: 






<gI< 6M? 
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(37) 



As we will see below, this inequality leads to a lower 
bound on the mass of the curvaton which noticeably 
shrinks the allowed region in the space of parameters. 

If we do not assume that inflation lasted for many more 
e-folds than those we have causal access to, then the limit 
above weakens. The inequality: 



< af < 6M? 
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(38) 



does not constrain the value of a.^, strongly and the al- 
lowed range for the mass of the curvaton is big. We are 
left with the single constraint m > T^ which this leads 
to a much smaller lower bound for the mass. We will 
see that this scenario is somewhat less interesting than 
the one in which m is more restricted because the for- 
mer calls for a certain amount of fine tuning if CDM is 
required not to leave any isocurvature trace. Thus, we 
treat this "short inflation scenario" as a caveat that we 
discuss at the end of this section, while we assume the 
long inflation scenario in what follows. 



B. Constraints coming from the power spectrum 

According to our assumptions, the curvaton must be 
responsible for most of the power at the pivot scale. Us- 
ing (111 and the relation of r to /nl, we can find an 



additional constraint on the parameters that determine 
the density contrast. Equation (111 becomes: 
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Solving for ct* we get: 



-ff* 



127r/NL 4.8 X 10-5 



= 27.6iJ, 
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(39) 



(40) 



The lower bound on cr* , ( 36 1 , provides a powerful relation 



between the energy scale for infiation and the mass of the 
curvaton: 



7 X 10-Vnl < tt < 1 



(41) 



where the upper bound comes from the requirement that 
the curvaton starts oscillating after the pivot scale ex- 
its the horizon. As a matter of fact, this bound is even 
tighter because the curvaton should not oscillate before 
the end of inflation, and Hend < H^. However, the 
precise relation between these two values of the Hubble 
factor depends on the particular inflationary potential. 
Since we wish to keep the analysis valid for a generic po- 
tential, we simply use (41 1 keeping in mind the previous 
statement. 

We also considered bounds coming from the tilt of the 
spectrum. Even though it is not explicitly quoted in [2], 
looking at the two dimensional probability distribution of 
a and Us in Fig. 9 of that paper, we estimate a marginal- 
ized 2 sigma region of: 0.93 < n^ ^ 1. The spectral index 
in this model is related to the primordial parameters by 
i4j: 



n,, = 1 






(42) 



We must make sure that the two additional factors in 
the RHS of this equation are not in conflict with the 
experimental bound. 

The second factor, is directly related to the ratio of 
primordial tensor to scalar amplitude [23]. Due to the 
low energy during inflation, the curvaton model predicts 
non observable primordial tensor modes, which implies 
T < 0.07 [2^. Thus, inside this workframe, it is expected 

that the ratio is limited from above by 2jj^ < 0.01 which 
would be irrelevant given the current precision in the 
determination of the spectral tilt. 

On the other hand, in order to prevent the generation 
of a blue tilt in conflict with observations, we impose the 
mild requirement that: 



m 



^il <''<"'"'> H, 



<0.1 



(43) 



which reduces the allowed range \41\ (this was also 
pointed out in [J]). 

Plugging in our latest relations into the decay rate, we 
get: 



r. = 6 • m 



1.59 X W 



(44) 



where 6 is a number that parametrizes the curvaton mass, 
m = b ■ H^,, and can take any value inside the interval 
(7x10-5/nl,0.1). 

After all the considerations, we can see that T^ de- 
pends mainly on the scale of inflation, and only mildly 
on the precise value of the mass of the curvaton. Indeed, 
the corresponding decay temperature, only depends on 

Vb. 



C. Mass of the CDM candidate particle 



We stick to the convention described in section III CI 
regarding the definition of the time of CDM creation. For 
candidates such as WIMPS, their creation corresponds 
to their freezing-out of thermal equilibrium in the early 
universe. In this case, we use the relation 1251: 



t cdm 



Wcdn 

20 



"^cdr 



<20v/Tviij 
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(45) 



to find an upper bound to the mass of the CDM candi- 
date. Note that the mass of the heaviest thermal relic 
corresponds to rn-cdm — 240 TeV [55]. Heavier dark mat- 
ter species must have been created out of equilibrium, in 
order to not over close the universe. In this case, the rela- 



tion of the mass to the creation temperature is not (45 1, 



and the limits apply to the temperature of the universe 
at the era of their creation. 

For simplicity, during this section we will assume that 
the non-gaussian contribution to the primordial power 
spectrum is /nl — 100. The results are easily generaliz- 
able to a different value of /nl, as long as it fulfills the 



large-non-gaussianity requirement specified in section III 



Also, we show below that the results are pretty robust 
against variations of this parameter. In order to mini- 
mize the size of the space of parameters, we fix the value 
of & = 0.1 and check at the end that the effect of this 
factor on the results is negligible. Note that in any case, 
this fixed value of b will generate the most conservative 
results as it sets the highest possible upper bound for Fo-. 
In order to prevent a neutrino isocurvature mode [9], 
the decay of the curvaton must occur before the neutrino 
decoupling era, at T^, ~ IMeV: 

F^ > — ^ ^ iJ* > 1.4 X lO^GeV (46) 

which leaves a narrow margin for the energy scale of in- 
fiation: 

6.7 X lO^^GeV > Y^^'^ > 2.4 x lO^GeV (47) 

as was pointed out in [27j (note that this range could be 
shrunk further if e ^ 1, but this depends on the partic- 
ular inflationary potential) . Combining equations ( 45 ) 



and ( 44 ) , we arrive at a new relation that connects infla- 



tionary phenomenology and particle physics beyond the 
standard model: 



TOcdm <3.2x lO^AfpL^ii^; 



(48) 



It applies for masses smaller than ^ 10^ GeV, that is, for 
H^ < 9.6 x 10^° GeV. For higher values of the Hubble 
factor, the bound is imposed on Tcdm: 



Tcdm <1.6x lO^Mp/ff. 



(49) 



with an absolute maximum for Tcdm is T'cdm < 1.8 x 10^ 
GeV. 



In Fig. [T] we plot these bounds as a function of the 
inflationary scale, i?*. The shaded areas above the solid 
line correspond to values of the parameters that would 
generate an isocurvature mode, for a curvaton mass of 
TO = 0.1 X H^,. The dotted line limits the region in which 
the creation temperature of CDM can be directly related 
to its mass via equation ( 45 1 . Above this line we can only 



constrain the temperature of the universe at the epoch 
of CDM creation. To assess the generality of the model, 
we also plot the limit resulting for the minimum possible 
mass of the inflaton, corresponding to 6 = 7 x 10"'^ and 
see that the conclusions are practically the same (dashed 
line). 
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FIG. 1: Upper bound on the mass of a CDM candidate as a 
function of the inflationary scale. We used the two bounding 
values for the mass of the curvaton: m'""^ = O.IH* (dashed), 
and m"""" = 7 • 10"^i/* (solid). The dotted line corresponds 
to TTicdm = 2.4 x lO^GeV, the value above which the bounds 
apply to the temperature of DM creation. The vertical dot- 
dashed line shows the tighter constraint on the scale of in- 
flation coming from assuming a quadratic potential for the 
inflaton. The value of /nl has been fixed to 100. 



Note that the consequences of the eventual detection 
of a CDM particle of a mass of around 100 GeV, would 
contract the allowed range by two orders of magnitude, 
leaving a narrow range of ~ 4- 10^ — 10^^ GeV for the pos- 
sible energy scale for inflation. It is also remarkable that 
once one assumes a particular shape for the inflaton po- 
tential, the parameter e can be determined. Then, since 
its value is generally smaller than 1, the upper bound 
for the scale of inflation will lower, shrinking further the 
allowed region. In particular, for a quadratic potential, 
we would flnd e ~ 10~^, and if 50 (or the maximum value 
for 7J*) would be one order of magnitude smaller. We 
show this with a dot-dashed vertical line in Fig. [l] . In 
this case, the absolute upper limit for the creation tem- 
perature of a CDM candidate is, Tcdm < 1-7 x lO*" GeV. 

The two temperature upper bounds would be relevant 
for primordial black hole (PHB) production, although 
the evolution and flnal density of such a species is highly 



dependent on the particular details of each model. Fur- 
thermore, small size, quickly evaporating PBH's could be 
generated at very early times leading to a large entropy 
release (see, e. g. [28]) that would modify the conclu- 
sions of this analysis. Therefore the possibility of PBH 
generation and the compatibility with the model, is left 
for future work. 

Gravitationally generated wimpzillas are created 
mostly at the end of inflation, for inflationary scales of 
Hend ~ 10^^ GeV [29j. This corresponds to a tempera- 
ture above the upper bound, and thus they are incompat- 
ible with this curvaton model. Nevertheless, wimpzillas 
can also be generated during reheating and preheating 
leading to model dependent masses or temperatures of 
creation. The existence of these second kind of wimpzil- 
las would not be in conflict, in general, with the curvaton 
bounds. 

"Short inflation" scenario: If inequalities (37 1 no 



longer apply, the lower bound on m comes from demand- 
ing: 



r„ 



m • H^ 



1.59 X 10^ 



> T'^/M, 



PL 



(50) 



which becomes the most extreme for the maximum value 
of H^,, in eq. (l7|. In that case, the above inequality is 
fulfilled as long as: 



m > 1.14 X 10^2^ AfpL 



(51) 



which corresponds to a value of 6 ~ 2 x 10~^^. Essentialy, 
extending the allowed range for b opens up the space of 
parameters and the relation between Tcdm Sind the scale 
of inflation is incomplete without the speciflcation of the 
value of b. However, it is significant that the range opens 
up towards its lower end since allowing b to take smaller 
values only narrows down the allowed region of Fig. IT] 
This leaves less and less time for the CDM to freeze out 
after the decay has occurred and would lead ultimately 
to a CDM that has to decouple just before neutrinos do. 
This makes the model fine tuned as, to some extend, less 
likely because the masses of the curvaton and the CDM 
particle which are unrelated a priori, must conspire to 
prevent the generation of isocurvature modes. 

Now we study how the results change when we drop 
the assumption that /nl = 100. In Fig. [2] we plot the 
upper bound on rricdm as a function of /nl for several 
fixed values of the energy scale. The range used for the 
possible value of the non-gaussianity is: 

10 < /nl < 400 

where the upper bound is an approximation of the result 
in [37], where they set the limit /nl < 522 • r^/^ for the 
curvaton model of inflation. We see how, for a fixed infia- 
tionary scale, the maximum allowed Tcdm hardly varies 
with /nl- This is again due to the strong dependence of 
^cdm on the scale of inflation, while it only depends lin- 
early on the non-gaussianity factor. With this we show 
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FIG. 2: Upper bound on the mass of a CDM candidate as 
a function of non-gaussian fraction. We used three different 
fixed values for the scale of inflation: 10^ GeV (solid), 10^" 
GeV (dashed), 10^^ GeV (dotted). The value of b used is 
6 = 0.1 



that the analysis is robust against variation of the non- 
gaussianity contribution. This is particularly convenient, 
given the existing uncertainties in the estimation of /nl- 



V. CONCLUSIONS 

Given the recent advances on the detection of primor- 
dial non-gaussianity in the CMB, we examined the impli- 
cations of an eventual observation of a large non-gaussian 
signal (/nl ~ 100) for a very simple curvaton model. 
In particular, we studied the residual isocurvature signal 
that could be imprinted in this scenario. 

It is found that, if the curvaton is indeed the mecha- 
nism responsible for seeding a primordial spectrum with 
an important non-gaussian contribution, it cannot induce 
as well any trace of primordial isocurvature. This leads to 
two important conclusions that help confining the model: 
1. The curvaton field cannot decay into dark matter, 2. 
The dark matter species must decouple or be created af- 
ter the curvaton has decayed. 

Conclusion 2, is particularly interesting from a particle 
physics point of view: 

Firstly, it allows us to set an absolute upper bound for 
the creation temperature of dark matter, which has to 
be lower than 1.7 x 10^ GeV for the case where the infia- 
tionary potential is quadratic. Gravitationally produced 
wimpzillas are in clear conflict with this bound and the 
assumptions of the model. However, due to the strong 
model dependence of non-thermal dark matter genera- 
tion, it is not straightforward to rule out any other class 
of such species. 

Secondly, it implies a relation between the temperature 
of creation of a cold dark matter candidate and the scale 



of inflation: 



Tcdm< 1.6 X W^Mpr^m 



V/nl 



(52) 



which hnks the physics of the early universe to the physics 
of the dark matter sector. This relation is shown to be ro- 
bust against variations of the non-gaussian component, 
and would be specially significant if a WIMP was ob- 
served in future experiments, since a portion of the space 
of parameters for this model would be ruled out. 

We remark the importance of the conclusions drawn 
above, because if the curvaton model turns out to be 
relevant inside the inflationary picture, this connection 



to dark matter physics could be the only probe into the 
scale of inflation. 
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